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ABSTRACT 



Aims. The interactive software package iCosmo, designed to perform cosmological calculations is described. 

Methods. iCosmo is a software package to perform interactive cosmological calculations for the low redshift universe. The com- 
putation of distance measures, the matter power spectrum, and the growth factor is supported for any values of the cosmological 
parameters. It also performs the computation of observables for several cosmological probes such as weak gravitational lensing, 
baryon acoustic oscillations and supemovae. The associated errors for these observables can be derived for customised surveys, or 
for pre-set values corresponding to current or planned instruments. The code also allows for the calculation of cosmological forecasts 
with Fisher matrices which can be manipulated to combine different surveys and cosmological probes. The code is written in the 
IDL language and thus benefits from the convenient interactive features and scientific library available in this language. iCosmo can 
also be used as an engine to perform cosmological calculations in batch mode, and forms a convenient evolutive platform for the 
development of further cosmological modules. With its extensive documentation, it may also serve as a useful resource for teaching 
and for newcomers in the field of cosmology. 

Results. The iCosmo package is described with various examples and command sequences. The code is freely available with docu- 
mentation at http : //www . icosmo . org along with an interactive web interface and is part of the Initiative for Cosmology, a common 
archive for cosmological resources. 
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1. Introduction 

Cosmology is a rapidly advancing field thanks to recent progress 
in instrumentation, numerical simulations and theoretical meth- 
ods. The last few decades have provided a previously inacces- 
sible wealth of data with the advent of large surveys of the 
cosmic microwave background, galaxy clust ering, weak gravita- 
tional lensing, cluste rs and supernovae (see ISpergel et al.ll2007t 

. iDunklev et al.ll2008l for recent results and references therein). 

■ One of the challenges of cosmology is to link the observables 
which can be derived from these cosmological surveys with 
the parameters of an underlying cosmological model. Another 
necessary task is to predict the performance of future surveys 
and future observational probes in order to maximise the po- 
tential of planned experiments. Existing publicly available com- 
putational tools (Seljak & Zaldarriaga 1996; Lewis et al. 2000t 
iLewis & Bridld 120021) are widely used within the cosmologi- 
cal community, yet there does not exist at present a unified 
tool which considers the low redshift universe and its associated 
probes. 

In this paper, we present iCosmo, a software package to per- 
form interactive cosmological calculations. The code can also 
be used as an engine to perform predictions for any cosmological 
model and also forms a convenient platform for the development 
of further cosmological modules. With its extensive documenta- 
tion, it may also serve as a useful resource for teaching and for 
newcomers in the field of cosmology. 

The code performs the calculation of distance measures, the 
linear and non-linear matter power spectrum, the growth factor, 
volume elements and other quantities for any CDM cosmolog- 
ical model. It also performs the computation of observables for 



several cosmological probes such as weak gravitational lensing, 
baryonic acoustic oscillations (BAO) and supernovae. The asso- 
ciated eiTors for these observables can be described for arbitrary 
survey parameters, or alternatively for a number of pre-set val- 
ues corresponding to planned surveys and experiments. The code 
also provides forecasts for future cosmological surveys, by al- 
lowing the user to compute, manipulate and combine Fisher ma- 
trices. The code is written in the IDL language and thus benefits 
from the simple syntax, large scientific libraries and convenient 
interactive plotting environment of this language. In particular 
all the above calculations can be completed in only a few lines 
of code. 



The code is freely available at 'http : //www . icosmo . org' 
with full documentation and help files. It is part of the Initiative 
for Cosmology website, a global archive for cosmological re- 
sources. An interactive web interf ace to iCosmo is also av ailable 
at this website (see description in lKitching & et aLll2008l) . 

In this paper, we describe the main features of iCosmo and 
provide a number of examples illustrating its use. In section |2] 
we describe the main conventions and architecture of the code. 
In section [3] we describe how the code can quickly calculate the 
evolution of the basic cosmological quantities, such as distance 
measures and growth factor, as a function of redshift. In sec- 
tion |4] we describe how iCosmo calculates the observables for 
several cosmological probes such as weak lensing, BAO and su- 
pernovae. Section |5] shows how the code can compute and ma- 
nipulate Fisher matrices. Our conclusions and possible future de- 
velopments are presented in section|6l 
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2. General presentation 

The iCosmo code is a package of routines written in the IDL 
language. It is divided into several directories corresponding to 
the different levels of the calculations: 

- general , plotting : general utility and plotting routines. 

- cosmo, expt: routines to define fiducial cosmological and 
survey parameters and compute basic cosmological quanti- 
ties. 

- lensing, bao, sne: routines to compute weak lensing, BAO 
and supernovae observables. 

- fisher: routines to create and manipulate Fisher matrices 
for cosmological forecasts. 

The routines in each of the directories make use of the vari- 
able structures in IDL. Their functionality is indicated by the 
first few letters of their name with the following conventions: 

- set: set up basic structures with calculation variables 

- mk: make or compute a new structure using the information 
in an input structure 

- get: get a substructure or derived quantity from an existing 
structure 

- pit: plot 

- rd, wt: read and write data from and into a file 

A description of the main iCosmo routines and associated 
directories can be found in table [1] Detailed explanations about 
any of the routines can be obtained by typing icosmojielp , 
'routine name' at the IDL command line. A readme file is avail- 
able in the distribution and gives installation instructions and a 
quick start tutorial. 

Table |2] shows an example of a typical call sequence for 
iCosmo, with the different levels of the code separated with hor- 
izontal lines. The following sections describes the use of each 
levels of the code using examples of call sequences and output 
figures. 

3. Basic cosmological quantities 

The first level in iCosmo defines the cosmolo gical model and 
computes basic cosmological quantities (see e. g iPeacockl i 1 999t) 
for their definitions, and Refregier et al. (2004) for conventions). 
The first step is to define the fiducial parameter structure using 
the set_f iducial function (see Inst. 1 in Table|2]). The resulting 
fid structure contains the following three substructures, which 
contain the following information: 

- cosmo: input cosmological parameters 

- expt: experiment parameters, such as survey area, number 
of galaxies or SNe as a function of redshift, or pre-set sur- 
veys. 

- calc: calculation parameters, such as redshift and ^-range, 
fitting functions, etc. 

These parameters have default values which can easily be mod- 
ified using the syntax of instruction 1 in Table |2] 

The next step is to compute all the basic cosmological quan- 
tities using the mk_cosmo function (see Inst. 2 in Table |2]l. The 
resulting cosmo structure is organised into several substructures: 

- const: all original fiducial quantities as well as constant 
quantities derived from the fiducial parameters (e.g. the 
Hubble radius) 

- evol: Evolving scalar quantities such as distance measures, 
the growth factor, etc, tabulated as a function of redshift 



- pk: Linear and non-linear matter power spectrum P(k, z) tab- 
ulated as a function of wavenumber k and redshift z 

The evolving scalar quantities can then be easily plotted us- 
ing the plt_cosmo routine (see Inst. 3 in Table|2]i. As an exam- 
ple, the following call sequence 

) fid=set_f iducial (cosmo={w0 : -Q . 9) , calc={f it_tk : 1)) 

) cosmo=mk_cosmo(fid) 

) plt_cosmo,cosmo, 'z' , 'da' 

produces Figure [T] which shows the angular diameter distance 
Da(z) as a function of redshift z- The fiducial cosmology is 
(h = 0.7, Qi = 0.045, n,„ = 0.3, Qa = 0.7, wq = -0.95, w„ = 
0.0, n - 1.0, T = 0.09, cTg = 0.8), unless otherwise stipulated. 
For example, the call sequence above changes wq to -0.9, which 
is used in the rest of this paper 

The plotting routine plt_pk can be used to plot the power 
spectrum P{k, z) as a function of wavenumber A: at a redshift z- 
For example. Figure |2] shows the linear and non-linear power 
spectrum at z = and 1 and is produced by the additional 
sequence 

) plt_pk, cosmo, z=0,xran= [0.001, 10.] 

) plt_pk , cosmo , z=0 , /over , /linear , linestyle=2 

) plt_pk, cosmo, z=l, /over, color=2 

) plt_pk , cosmo , z=l , /over , /linear , linestyle=2 , color=2 
Note that the keyword over can be used to overlay several 
results on a single figure. 
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Fig. 1. Example of an iCosmo output showing the angular- 
diameter distance Da(z) as a function of redshift. The plot was 
derived using the plt_cosmo routine. See text for command se- 
quence. 



4. Cosmological observables 

The next level of iCosmo consists of the calculations of the ob- 
servables related to different cosmological probes. At present, 
weak lensing, baryon acoustic oscillations and supernovae 
are supported. The former probe is imple mented usirig shea r 
power spectrum tomography as described in Hu & J ainI (l2004l) : 
Refregier et al. (2004); Amara & Refregier (2007). For BAO, ra- 
dial and tange ntial distance scales , along with their errors are de- 
rived from the IB lake et al.l (|2006|) fitting formulae. We calculate 
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Directory 


Routine 


Description 


cosmo 


set_fiducial 


Create fiducial parameter structure 




mk.cosmo 


Compute basic cosmological quantities 




get_pk 


Extract 3D power spectrum 


expt 


mk.survey 


Compute survey parameters 


lensing 


mk_cl_tomo 


Compute lensing power spectrum 




mk_cl_cov_tomo 


Compute covariance errors 




mk_fisher_lens 


Compute weak lensing Fisher matrix and errors 


bao 


mk.bao 


Compute BAO distance measures 




mk_bao_cov 


Compute BAO covariance errors 




mk_fi sher.bao 


Compute BAO Fisher matrix and errors 


sne 


mk.sne 


Compute SNe magnitude-redshiit relation 




mk_sne_cov 


Compute SNe covariance errors 




mk_fisher_sne 


Compute Sne Fisher matrix 


fisher 


checkjnatrix 


Check matrix is positive definite 




comb.fisher 


Combine two Fisher matrices 




margin.fisher 


Marginalise over unwanted parameters 


plotting 


plt_cosmo 


Plot cosmology parameters 




plt_pk 


Plot 3D matter power spectrum 




plt_sv 


Plot survey properties 




plt_cl 


Plot weak lensing correlation function 




plt_bao 


Plot BAO distance measures 




plt_sne 


Plot SNe magnitude-redshift relation 




plt_fisher 


Plot Fisher matrix errors (combined ID and 2D errors) 




plt_fisher_lp 


Plot ID likelihood errors 




plt_fisher_2p 


Plot 2D error ellipses 



Table 2. Example of a typical call sequence for iCosmo. The different levels of the code are separated by horizontal lines. 





Instruction 


Description 


1 


f id=set_f iducial (cosmo_in={omega_m : . 3 Ij) 


Set fiducial parameters, choosing a non-default value of the parameter Q.,„ 


2 


cosmo=mk_cosmo(fid) 


Compute basic cosmological quantities (distance measures, power spectrum, etc) 


3 


plt_cosmo , cosmo , ' z ' , ' da ' 


Plot angular-diameter distance Z)a(z) as a function of redshift z 


4 


plt_pk , cosmo , z=8 


Plot linear matter power spectrum P(K) at z = 


5 


sv=mk_survey(fid, 'svl') 


Read survey parameters for the pre-set DUNE experiment 


6 


cl=mk_cl_tomo(f id , cosmo , sv) 


Compute weak lensing power spectrum C; 


7 


cl_cov=mk_cl_cov_tomo(f id , cl , sv) 


Compute weak lensing power spectrum errors 


8 


plt_cl , cl , [8,8] , cl_cov , / errors 


Plot lensing power spectrum and associated errors 


9 


fish=mk_fisher_lens(fid, sv) 


Compute fisher matrix for the specified lensing survey 


10 


margin_f isher , f ish, f ish2 , [0 , 1 , 1 , , 8 , 8 , 8 , 0] 


Marginalise fisher matrix over unwanted parameters 


11 


plt_fisher_2p, fisher2 


plot fisher matrix error constraints for interesting parameters wo and 



magnitude -redshift relation for s u pernovae using the forma lism 
outlined in tegm ark et all ( ll998l) : lHuterer & TurneJ (1200 Ih and 
include both intrinsic scatter and observed variance. 

A typical call sequence is given in Instructions 5-8 in Table|2] 
The sequence first defines the survey using mk_survey and then 
computes the observable using mk_cl_tomo, mk_bao, mk_sne, 
for lensing tomography, BAO and SNe respectively. The statisti- 
cal errors in these observables for the fiducial survey can be de- 
rived using the corresponding mk_cov routines listed in table [T] 

As an example, the following prod uces Figure [3] a plot of 
weak lensing power spectra for DUNE (lRefregierll2008h with its 
associated Icr error bars: 

) fid=set_fiducial('DUNE' , exptj.n={svl_n_zbin: 2)) 

) cosmo=mk_cosmo(£id) 

) sv=mk_survey(£id, ' svl ' ) 

) cl=mk_cl_tomo (fid, cosmo, sv) 



) cl_cov=mk_cl_cov_tomo(f id , cl , sv) 

> plt_cl,cl, [1, 1] ,cl_cov, /errors, yran=[le-7, le-3] 

) plt_cl,cl, [0, 1] ,cl_cov, /errors, /over, linestyle=l 

) plt_cl,cl, [0,0] ,cl_cov, /errors, /over, linestyle=2 

The different power spectra correspond to two different to- 
mographic redshift bins with median redshifts 0.68 and 1.36, as 
well as their cross-power spectra. Similar call sequences can pro- 
duce predictions and errors for the BAO distance measures and 
for SNe Hubble diagrams. 



5. Fisher Matrices 

The final level of iCosmo is related to the computation and ma- 
nipulation of Fisher matrices to assess the constraints on cos- 
mological parameters that can be achieved with future surveys. 
For this purpose, fisher matrices can first be computed using the 
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Fig. 2. Matter power spectrum at z = and 1 , for top and bot- 
tom lines respectively. Linear and non-linear power spectra are 
shown as dashed and solid lines, respectively. The plot was de- 
rived using the plt_pk routine. 
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Fig. 3. The weak lensing power spectrum and associated Icr er- 
ror bars for the DUNE weak lensing survey for the first (upper 
line) and second (middle line) tomographic bins (with median 
redshifts of 0.68 and 1.36 respectively) as well as the cross- 
correlation between the bins (lower line). 



mk_£isher routines for each probe (see Table[T]i. The computa- 
tion is done by taking excursions from the fiducial cosmologi- 
cal model for each of the parameters of interest. Fisher matrices 
can be combined using the combine_fisher routine to derive 
joint cosmological constrains with several surveys and probes. 
The routine margin_fisher can be used to marginalise over 
unwanted parameters. The resulting constraints can be plotted 
using the plt_fisher routines. 



A typical call sequence to compute and plot fisher matrices 
can be found in instructions 9-11 in Table|2] As an example, the 
following sequence computes the fisher matrices for a DUNE 
weak lensing survey, a half-sky BAO spectroscopic survey 
(20000 deg^), and their combination. 

) fid_lens=set_£iducial('DUNE') 

) £id_bao=set_£iducial( 'bao±alfsky ' ) 

) sv_lens=mk_survey(£id_lens, 'svl') 

) sv_bao=mk_survey (£id_bao , ' sv2 ' ) 

) f_lens=mk_fisher_lens(fid_lens , sv.lens) 

) f_bao=mk_£isher_bao(fid_bao , sv_bao) 

) f_comb=comb_£isher (£_lens , £_bao) 

The following call sequence produces figure |4] showing the 
68%CL constraints on the dark energy parameters Q.de and wo 
expected for each surveys separately (blue and red ellipse) and 
jointly (solid green ellipse). 

) margin_£isher , £_lens , f_lens2 , [Q , 1 ,Q ,6) ,0,0,8 , 1] 

> margin_£isher ,£_bao, f_bao2 ,[0,1,0,0,0,1] 

) margin_£isher , £_comb, £_comb2 ,[0,1,0,0,0,0,0,1] 

) plt_£isher_2p , £_bao2 , /no£ill , color=2 

) plt_£isher_2p , £_lens2 , /no£ill , linestyle=2 , /over 

) plt_£isher_2p , £_comb2 , /over , color=3 

The margin_fisher routine was used to retain only the pa- 
rameters to be plotted, and marginalise over all the other param- 
eters. 
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Fig. 4. Marginalised 68%CL constraints on the dark energy pa- 
rameters Q.DE and wq expected for the DUNE weak lensing 
(blue), a full sky BAO survey (red) and their combination (solid 
green). This figure was derived using the Fisher matrix routines 
of iCosmo. 



6. Conclusions 

The iCosmo package provides a convenient and flexible inter- 
active tool to make cosmological calculations and can also be 
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used as a computation engine in batch mode. Using all the built 
in IDL libraries, it is also a convenient platform for the de- 
velopment of further cosmological routines a nd for teaching. 
The c ode, along with an interactive web tool ( Kitching & et alj 
12008) and various cosmological resources are freely available at 
pittp : //www . icosmo . org 

In the future, we plan to add a number of features to the 
code, such as the halo model, higher order clustering statistics, 
and interfaces with CMB Boltzmann codes. Contributions from 
the community are encouraged. 

Acknowledgements. The authors thank the many people who have contributed to 
iCosmo over the years, pailicularly, Richard Massey, David Bacon, Joel Berge, 
Florian Pacaud, Romain Teyssier, Jean-Baptiste Juin, Dominique Ivon, and Ivan 
Debono. 



References 

Amara, A. & Refregier, A. 2007, MNRAS, 381, 1018 
Blake, C, Parkinson, D., Bassett, B., et al. 2006, MNRAS, 365, 255 
Dunkley, J., Komatsu, E., Nolta, M. R., Spergel, D. N., & et al. 2008, ArXiv 
e-prints 

Hu, W. & Jain, B. 2004, Phys. Rev. D, 70, 043009 
Huterer, D. & Turner, M. S. 2001, Phys. Rev. D, 64, 123527 
Kitching, T. & et al. 2008 

Lewis, A. & Bridle, S. 2002, Phys. Rev., D66, 103511 

Lewis, A., Chalhnor, A., & Lasenby, A. 2000, Astrophys. J., 538, 473 

Peacock, J. A. 1999, Cosmological Physics (Cosmological Physics, by John 

A. Peacock, pp. 704. ISBN 052141072X. Cambridge. UK: Cambridge 

University Press, January 1999.) 
Refregier, A. 2008, Experimental Astronomy, 26 
Refregier, A., Massey, R., Rhodes, J., et al. 2004, AJ, 127, 3102 
Seljak, U. & Zaldaniaga, M. 1996, Astrophys. J., 469, 437 
Spergel, D. N., Bean, R., Dore, O., Nolta, M. R., & et al. 2007, Apjs, 170, 377 
Tegmark, M., Eisenstein, D. J.. & Hu, W. 1998, ArXiv Astrophysics e-prints 



